Coherence scale of the two-dimensional Kondo Lattice model. 
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A doped hole in the two-dimensional half-filled Kondo lattice model with exchange J and hopping 
t has momentum (71", n) irrespective of the coupling J/t. The quasiparticle residue of the doped hole, 
tracks the Kondo scale, Tk , of the corresponding single impurity model. Those results stem 
from high precision quantum Monte Carlo simulations on lattices up to 12 x 12. Accounting for 
small dopings away from half-filling within a rigid band approximation, this result implies that 
the effective mass of the charge carriers at the Fermi level tracks 1/Tk or equivalently that the 
coherence temperature T co h oc Tk- This results is consistent with the large-N saddle point of the 
SU(N) symmetric Kondo lattice model. 

PACS numbers: 71.27. +a, 71.10.-w, 71.10.Fd 



In a Fermi liquid, the coherence scale corresponds to 
the energy scale below which thermodynamic properties 
are determined by the Fermi surface: the specific heat is 
linear in temperature and the spin charge uniform suscep- 
tibilities temperature independent. This scale is inversely 
proportional to the effective mass of the charge carriers. 
In this article, our aim is to compute the coherence scale 
in models of heavy fermion materials 0. Our starting 
point is the Kondo Lattice model (KLM), 
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on a square lattice. The first term describes a band 
of non-interacting electrons in a tight binding approxi- 
mation (cl creates an electron in a Wannier state cen- 

tered around lattice site i with z-component of spin a) 
and the sum runs over nearest neighbors. The spin 
degrees of freedom of the conduction electrons, S° = 

^y~! 00 'ct a s s'd , with a the Pauli spin matrices, in- 

^ £ — ' b > b 2,S ' i,s 

teract antiferromagnetically, J > 0, with a lattice of spin 
1/2 magnetic impurities, St. The KLM at J/t « 1 
stems from the periodic Anderson model (PAM) in the 
limit where charge fluctuations are negligible. 

In the limit of a single impurity, the model is well un- 
derstood: at high temperatures the impurity spin is es- 
sentially free and at low temperatures it is screened via 
the formation of a many body singlet state of the con- 
duction electrons and impurity spin. The characteris- 
tic, universal, energy scale describing this crossover from 
the free to the screened magnetic impurity is the Kondo 
temperature, Tk 0. I n the lattice case, it is appeal- 
ing to view the heavy fermion metallic state as a conse- 
quence of a coherent, Bloch-like, superposition of indi- 
vidual Kondo screening clouds. The coherence scale of 
the metallic state has been investigated in details within 
the large-N approximation for the KLM Q and the dy- 
namical mean field approximation for the related PAM 
|j| ■ Both approaches yield T co h cx Tk as a function of J / 1 
with a proportionality factor depending strongly on the 
band-filling. Since both approximations neglect spatial 
fluctuations, they do not capture the Ruderman-Kittel- 



Kasuya-Yosida (RKKY) interaction which introduces a 
new energy scale into the problem and drives the ground 
state through a magnetic order-disorder quantum phase 
transition. The new result in this Letter is that the re- 
lation T co h cx Tk still holds when the two-dimensional 
KLM is solved numerically exactly by means of quantum 
Monte Carlo (QMC) simulations. 

To tackle this problem, we have to adopt an indirect 
route since the sign problem inhibits simulations away 
from half-filling. We hence investigate the problem of a 
single-hole doped into the Kondo insulating state at half- 
filling and then assume a rigid band to deduce the prop- 
erties of the metallic state at weak dopings. In the follow- 
ing, we describe our approach for the SU (N) symmetric 
model where in the limit N — > 00 a paramagnetic saddle 
point approximation becomes exact. We then compare 
the obtained results with QMC simulations for the SU (2) 
model of Eq. ^ 

Mean Field. The mean field we consider neglects 
magnetic fluctuations triggered by the RKKY interac- 
tion but describes the Kondo effect as well as the for- 
mation of the heavy electron state. The decoupling is 
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based on the equation: SSSi 

Di = E CT iju- Here, M = \ £ s>s , fU s , s > and 
/J are fermionic operators. At the mean-field level one 

adopts a real order parameter, r — (D^), and the con- 
straint of a single charge per f-site is imposed on aver- 
age via a Lagrange multiplier A. This mean field theory 
has been discussed extensively in Here, we will con- 
centrate on the J dependence of scales for a half-filled 
particle-hole symmetric conduction band. This symme- 
try pins the chemical potential to [i = and A = 0. Single 
particle properties are derived from the Green function 
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with self-energy 
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and e(k) = —2t (cos(k x ) + cos(k y )). Assuming a flat den- 
sity of states, -D(e), of width W satisfying J deD(e) = 1 
and D(e) — D(—e) we can solve the saddle point equa- 
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tions to obtain: 



r(T = 0) 



W 



2Jsinh(W/2J)' 
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The order parameter r vanishes at the Kondo temper- 
ature Tk- At the mean- field level, this energy scale is 
independent on the impurity concentration and hence 
matches the Kondo temperature of the corresponding 
single impurity problem. With the above flat density 
of states we obtain: 

T K oc We~ w/J oc — [Jr(T = 0)f for J/t « 1. (4) 



persion relation Er = 



The functional form of the self-energy leads to the dis- 

e(k) - Je 2 (k) + (Jr) 



residue Z^ 
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hole away from half-filling. At T = and in 
weak-coupling limit, both the quasiparticle gap, A 
fi — -E^tt, and residue track the same scale: 
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A qp oc Z (7! ^ oc [Jr(T = 0)] 2 oc T K forJ/< « 1. (5) 

Under a rigid band assumption (see Discussion section) 
this results implies that the coherence temperature of the 
metallic state at small hole dopings away from half-filling 
tracks the the Kondo scale. At the mean-field level this 
statement may be checked explicitly by solving the mean- 
field equations in the metallic state Q. 

The mean-field theory has many caveats. The order 
parameter r does not vanish because the constraint is 
taken into account on average. Magnetism is which is en- 
ergetically favorable in the small J/t limit is not included. 
Hence it is a-priori not clear that the scales described by 
the mean-field approach survive when the model is solved 
exactly. 



-G 



0.1 



k = (it, it), — t , 8 x 8, (n) = 1 



J/t = 1.45 




FIG. 1: Zero temperature Green function denned in Eq. |S| 
Fitting the tail of this quantity to the form Z^e~ Aqv{k)T yields 
the quasiparticle residue and quasiparticle gap at wave vector 
k 

QMC. To obtain zero temperature properties, we solve 
the half-filled KLM with the projector auxiliary field 



QMC (PQMC) method. The application of the method 
to the KLM has been describe extensively in Ref. 0, 0| • 
The algorithm scales as (3N 3 where, in the framework 
of the PQMC, (5 is an projection parameter in imaginary 
time and N is the number of orbitals. Hence, large values 
of (3 may be used to guarantee convergence to the ground 
state, but on the other hand the N 3 scaling renders large 
lattice sizes difficult to access. From the imaginary time 
displaced Green function, we compute both the quasipar- 
ticle gap as well as quasiparticle residue: 
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Clearly high precision 
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since Z^ 

results are required to extract the residue. Fig. ^ plots 
typical raw data. To reach such precision we use a newly 
developed method for the calculation of G(k,r) within 
the PQMC 0. 

Numerical Results. In numerical simulations of the 
two-dimensional KLM, the RKKY interaction triggers 
a magnetic order-disorder quantum phase transition at 
J c /t — 1.45 ±0.05 0| . The question then arises if the an- 
tiferromagnetic ordering destroys Kondo screening and 
hence the possible appearance of a Kondo scale in the 
numerical data. We answer this question by analyzing 
the single particle spectral function. Assuming that well 
under J c /t the impurity spins are frozen into their antifer- 
romagnetic ordering, we can use the mean-field Ansatz, 
(Sl' z ) = im{e I( 5' 1 , which leads to a single hole disper- 



sion relation: E(k) = -Je 2 (k) + (J/4) 



This ap- 



proximation correctly reproduces the quasiparticle gap 
A gp oc J/4 in the small J limit (see Fig. |3J but fails to ac- 
count for the observed dispersion relation in the spectral 
function. This statement is based on comparison with 
QMC data: the thin vertical lines in Fig. [21 tracking the 



dispersion relation — y e 2 (fc) + (J/4) 2 , do not account for 

the low energy features around the (7r,7r) point. On the 
other hand, comparison with the mean-field approxima- 
tion of the previous section (bold vertical lines in Fig. 
121, allows the interpretation that the low- lying features 
in the vicinity of k = (tt, it) stem from Kondo screen- 
ing and are shifted to energy scales of order J/t due to 
the onset of magnetic ordering. We have checked that 
this functional form of the dispersion relation survives 
down to J/t = 0.2 on 8 x 8 lattices. In the magnetically 
ordered phase, a doped hole at k = (jr, it) can scatter 
off a gapless magnon with momentum Q — (7r, it) thus 
generating a shadow feature at k = (0,0). Upon close 
inspection, such features are seen in Fig. [21 The quasi- 
particle dispersion relation is at best described in terms 
of partial Kondo screening of the impurity spins. The 
remnant non-screened moment orders antiferromagneti- 
cally 

Our major concern is the quasiparticle residue at k = 
(-7T, 7r). In Fig. [31 we plot this quantity as a function of 
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FIG. 2: Single particle spectral function as obtained from an- 
alytically continuing the imaginary time QMC data with the 
Maximum Entropy method. The single particle occupation 
numbers, n? = (ct c r ) are listed on the left hand side 
of the plot, and correspond to the weight under each spec- 
trum: 7mg = j_ dujA(k, iS). The thin vertical lines track 
the dispersion relation of hole in a static staggered magnetic 

field (see text): -\Je 2 (k) + {J/4) 2 . The bold vertical lines 
correspond to the dispersion relation stemming from the self- 
energy of Eq. |21 with the order parameter r determined self- 
consistently. 



J/t for 8x8 and 12 x 12 lattices. At small values of 
J/t size effects become more and more important. For 
J/t = 0.4 and 0.8 we have extrapolated the data to infi- 
nite sizes by fitting the L = 4, 6, 8, 12 QMC results to the 
form a + b/L + c/L 2 . To compare this scale to the single 
impurity Kondo scale, we have carried out simulations of 
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FIG. 3: Scales as a function of J/t as determined from QMC 
simulations (see text). The data for the quasiparticle gap are 
a result of extrapolation to the thermodynamic limit. The 
solid line corresponds to the mean-field (see section Mean- 
field) value of the Kondo temperature. Within our resolution 
we cannot distinguish it from the QMC value. 



the Kondo model: 
H = 
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with a standard implementation of the Hirsch-Fye QMC 
algorithm In this algorithm the CPU time scales 

as (] 3 N a where N corresponds to the number of lattice 
sites and (3 to the inverse temperature. This scaling be- 
havior allows one to carry out simulations on arbitrarily 
large lattices but limits the accessible temperature range. 
Form the data collapse of the impurity spin susceptibil- 
ity (see Fig. 0} we can extract the Kondo temperature 
for our given band structure. Note that the so obtained 
Kondo temperature (triangles in Fig|3J) compares remark- 
ably well with the mean- field value (solid line in Fig. |3J). 
At small values of J/t (J/t = 0.4) the Kondo tempera- 
ture is too small to compute with the Hirsch-Fye algo- 
rithm and we have to rely on the mean-field result to 
compare with the lattice QMC results. 

The data of Fig. |3]show several features, 
i) The major observation is that within the considered 
range of coupling constants the QMC results are consis- 
tent with 



OC T, 



K ■ 



(8) 



ii) Due to the occurrence of antiferromagnetic order be- 
low J c /t = 1.45 ± 0.05, the quasiparticle gap tracks J. 
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FIG. 4: Impurity spin susceptibility, = 

dr{m{(r)m{(0)), for the Kondo model as determined 
with a standard implementation of the Hirsch-Fye algorithm. 
From the collapse of the data, we can estimate the Kondo 
temperature for the listed values of J/t. 



Even in the presence of this large quasiparticle gap, the 
Kondo-like features in the single-hole spectral function 
- flat band around the (7r, 7r) point with Zt^^-s oc Tk - 
survive. 

iii) Within our resolution and limitation in lattice size, we 
observe no anomaly in the quasiparticle residue across the 
magnetic order-disorder phase transition ( J c /t ~ 1.45). 

Discussion Under the assumption of a rigid band, the 
result of Eq. [S] has important implications. Since the 
quasiparticle gap is determined by the (±7r, ±7r) points 
in the Brillouin zone, infinitesimal hole doping away from 
half-filling yields a metallic state with Fermi surface con- 
sisting of hole pockets centered around the (±7r, ±7t) 
k— points. This leads to a Fermi surface with large Lut- 
tinger volume containing both conduction and impurity 
electrons. The effective mass of this Fermi liquid at in- 
finitesimal hole dopings is inversely proportional to the 
quasiparticle residue: 

to* oc — . (9) 

Z(n,Tr) 

The coherence temperature, T C oh, is inversely propor- 



tional to the effective mass. Hence Eq. [5] along with 
our QMC result of Eq. [8] is equivalent to 

T coh oc T K , (10) 

Let us comment on the implicit assumptions lying behind 
the above equation. 

i) Support for the rigid band assumption follows from the 
fact the Kondo insulator is adiabatically linked to band 
insulator realized in the non-interacting periodic Ander- 
son model. This stands in contrast to the Mott insulting 
state. 

ii) That the effective mass is inversely proportional to 
to quasiparticle residue implicitly assumes that the en- 
hancement of m* is driven by the frequency dependence 
of the self-energy. We can obtain support for this state- 
ment by computing the effective mass as the inverse 
curvature of the dispersion relation around k = (jr, it) 
thus capturing the mass enhancement origination from 
the frequency as well as momentum dependence of the 
self-energy. Within our accuracy, this quantity tracks 
l/Z(K t7V ) as a function of J/t thus confirming the validity 
of our assumption. 

In conclusion, we have carried out high precision QMC 
calculations of the half-filled KLM on a square lattice. 
We have shown that from weak to strong coupling and 
across the magnetic quantum phase transition, a doped 
hole has momentum (7r, 7r) and quasiparticle residue 
which tracks the Kondo scale of the corresponding sin- 
gle impurity model. Assuming a rigid band, this result 
leads to the conclusion that the J-dependence of coher- 
ence scale in the metallic state at small dopings away 
from half-filling tracks the Kondo scale. Since this result 
compares favorably with large-N and dynamical mean- 
field theories which omit spatial fluctuations, we can fol- 
low the idea that the coherence scale is insensitive to the 
RKKY magnetic scale. 
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